Three-dimensional magnetic flux-closure patterns in mesoscopic Fe islands 
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We have investigated three-dimensional magnetization structures in numerous mesoscopic 
Fe/Mo(110) islands by means of x-ray magnetic circular dichroism combined with photoemission 
electron microscopy (XMCD-PEEM). The particles are epitaxial islands with an elongated hexago- 
nal shape with length of up to 2.5 yLxn and thickness of up to 250 nm. The XMCD-PEEM studies 
reveal asymmetric magnetization distributions at the surface of these particles. Micromagnetic sim- 
ulations are in excellent agreement with the observed magnetic structures and provide information 
on the internal structure of the magnetization which is not accessible in the experiment. It is shown 
that the magnetization is influenced mostly by the particle size and thickness rather than by the 
details of its shape. Hence, these hexagonal samples can be regarded as model systems for the study 
of the magnetization in thick, mesoscopic ferromagnets. 
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I. INTRODUCTION 

The spatial confinement of the magnetization in a fer- 
romagnet of mesoscopic or nanometric size can have a 
dramatic impact on both the magnetic properties and 
the magnetization distribution j3i 13, S S 13 • Besides the 
technological importance of understanding and control- 
ling such finite-size effects in magnetic nanoparticles that 
can be used for magneto-electronic devices, the physics 
of the magnetization in confined structures is a beautiful 
and exciting research topic _54] . A precise theoretical de- 
scription of the magnetization distribution in particles of 
a size between about 100 nm and a few microns can be ob- 
tained in the framework of micromagnetism _6,, _7J . Gen- 
erally, the micromagnetic equations can only be solved 
numerically. Mesoscopic magnetic particles are too small 
for a meaningful description in terms of macroscopic ma- 
terial constants like the magnetic susceptibility. On the 
other hand, the simple macrospin model - which could 
be solved analytically Q - can only be applied for very 
small (^10 nm) magnetic particles 0- It loses its va- 
lidity in larger ferromagnets, where the magnetization 
distribution is strongly inhomogcncous. 

The subtle size-dependence of the magnetization [Tol 
[m Il2| is particularly important when the particle size 
is comparable to the so-called magnetic exchange length 
1^ ^3 J *• e- a material-dependent length scale that de- 
scribes the order of magnitude for the extension of mag- 
netic inhomogeneities like domain walls or magnetic vor- 
tices. Besides the size, also the particle shape is known 
to have a strong influence on the magnetic properties 
of nanostructures. Numerous studies on magnetic parti- 
cles of different geometries, like rings 0|, rectane 
nanowires and various other geometries 0, 1 
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have been reported in the last years. Apart from a few 
exceptions where three-dimensional magnetic structures 
have been investigated numerically {e.g. 0, 0, |20|'). 
experimentally [21| or both [2^, most studies on mag- 
netic nanostructures that have been published over the 
last years have focussed on particles in which at least 
one dimension is so small that the magnetization in the 
sample is either two-dimensional (thin-film elements) or 
one-dimensional (thin nanowires). Therefore, compared 
to thin film elements, where the influence of the par- 
ticle shape on the magnetic properties has been amply 
studied, not much is known about the effect of spatial 
confinement in three-dimensional magnetic particles. 

Simil arly, s everal types of magnetic domain walls 
[23I I24 . 125I [23 have been investigated thoroughly in ex- 
tended magnetic films |^. In these studies, the sample 
was essentially confined in one dimension (the film thick- 
ness). Recently, the influence of geometric confinement 
on domain walls has been studied in systems with two- 
dimensional |27| and one-dimensional 28] magnetization, 
and different ty pes of head-to-head domain walls have 
been predicted |29i | and observed jsof in thin magnetic 
strips or rings of different width. However, the effect 
of a mesoscopic spatial confinement of three-dimensional 
magnetic domain walls has hardly been investigated up 
to now. 

In this paper we report on three-dimensional magne- 
tization distributions in mesoscopic ferromagnets. As a 
model system, we have analyzed several monocrystalline 
Fe islands of different size (up to ca. 2.5 yum length 
and up to ca. 250 nm thickness). These particles are 
small enough to display pronounced finite-size effects 
of the magnetization, while they are sufficiently large 
and thick to sustain inhomogeneous, three-dimensional 
arrangements of the magnetization. The experimen- 
tal part of this study is performed by means of X-ray 
magnetic circular dichroism with photoemission electron 
microscopy (XMCD-PEEM). The direct comparison of 
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FIG. 1: (color online) (a) Atomic force microscopy of a Fe 
island, about Ifim long, (b) Schematic presentation of a typi- 
cal particle geometry. The angles a — 45° and (3 = 35.26° are 
fixed parameters that result from the monocrystalline struc- 
ture of the islands. The shape is uniquely defined by the 
points Ai, A2, Bi, B2 and the thickness t. Note that this set 
of parameters is subject to various constraints. For instance, 
the edges (j4i, ^2) and {Bi, B2) are parallel to each other (the 
x-direction) . 



high-resolution experimental observations with highly ac- 
curate micromagnetic computer simulations allows for an 
unambiguous interpretation of some complicated, unex- 
pected magnetization distributions. This analysis leads 
us to the conclusion that the Landau pattern, which is 
well-known from magnetic thin-film elements, is merely a 
simple variant of a general, more complex magnetization 
structure. 



II. SAMPLE PREPARATION 

Mesoscopic monocrystalline Fe islands have been 
grown on atomically flat Mo(llO) buffer layers of thick- 
ness 10 nm deposited on Al2O3(ll-20). The growth 
is performed with pulsed laser deposition (PLD) in ultra- 
high vacuum at 850 K. Under these conditions, Fe islands 
are formed by means of the Stranski-Krastanov growth 
mode. A cover layer of Mo [three atomic layers], fol- 
lowed by Al[3nm] has been added to protect the samples 
from oxidation. A more detailed description of the fab- 
rication of self-assembled Fe/Mo(110) islands is reported 
elsewhere Is^,!!^!- The islands have a hexagonal shape, 
a flat top and inclined lateral facets, all atomically flat, 
see Fig. ^). 

If the growth proceeded under thermodynamic equi- 
librium, the shape of all islands would be expected to be 
identical, and described by the Wulff-Kaischev theorem 
(see [s^ and references therein) . We indeed observe that 
the islands share a common set of features: the angle of 
inclination of the lateral facets, the horizontal top facet, 
as well as the direction of all edges. However, the shape of 
the particles is generally not symmetric and not uniform. 
The vertical and in-plane aspect ratios slightly differ from 
one to another. This may arise from kinetic effects, the 
influence of atomic steps on the Mo(llO) surface or resid- 
ual strain. The size ranges between about 500 nm to 2500 
nm in length and about 50 to 250 nm in thickness. Our 
observations suggest that the shape of the particles is 
well reproduced by taking into account only {001} and 



{110} facets. The examination of the constraints related 
to these facets leads to an idealized geometrical construc- 
tion scheme for the shape of these islands as described 
in Fig^ Using this scheme, each island is uniquely de- 
scribed by the four points Ai,A2,Bi,B2 and the thickness 
t. Approximate values for the thickness and the coordi- 
nates of these four points have been extracted from elec- 
tron microscopy images to construct models that were 
used as input for the numerical simulations. The den- 
sity of the islands on the substrate is low enough, so that 
generally speaking, the magnetostatic coupling between 
the samples is negligible. Obviously, a few exceptional 
cases can also be found, where two islands have formed 
in close vicinity. 

These self-assembled Fe islands can be regarded as a 
model system for the study of magnetization distribu- 
tions in three-dimensional, mesoscopic particles. Due to 
their high structural quality the Fe/Mo(110) islands are 
particularly suited for this investigation. A comparable 
quality would be diffult to achieve with lithographically 
fabricated samples, especially for particles of such ele- 
vated thickness. The well-defined shape and structure 
of the self-assembled islands ensures that the magnetiza- 
tion distribution in the particle is governed by the spatial 
confinement, and not by roughness effects or by other 
structural inhomogeneities. However, the consideration 
of these islands as model systems for the magnetization 
in thick, structured magnetic particles is only valid if the 
details of the particle shape (e.g. the inclination angle 
of the facets or the precise hexagonal shape) are not of 
decisive importance for the resulting magnetic structure. 
Else, each particle would have its own characteristic mag- 
netic structure depending on its shape, thus precluding 
an extraction of valuable information of general validity. 
As will be shown in section FVII BI computer simulations 
demonstrate that the aspect ratio, the size and the thick- 
ness are the most important parameters, while details of 
the shape play only a minor role concerning the overall 
magnetization structure. 



III. EXPERIMENTAL SETUP 

The experimental investigation of the magnetic struc- 
tures was performed with XMCD-PEEM 35). The prin- 
ciple of XMCD-PEEM is the spin- and helicity-dependent 
cross-section for photoelectron emission with circularly 
polarized X-rays. The lateral resolution of this tech- 
nique is now better than 30 nm '36] , so that detailed 
images of the domain structure in mesoscopic elements 
can be obtained. The experiments have been carried 
out at the Nanospectroscopy beamline of the ELETTRA 
synchrotron radiation facility ^7]- The combination of 
low-energy electron microscopy (LEEM) and PEEM at 
ELETTRA makes it possible to image both the mor- 
phology and the magnetic structure with the same in- 
strument. In the XMCD-PEEM mode, the samples have 
been illuminated with monochromatic, elliptically polar- 




ized X-rays in a micro-spot of 20 x 5 /zm^ (in the horizon- 
tal and vertical direction, respectively) . The local photo- 
electron current density was imaged with PEEM. The 
photoelectron density is proportional to the magneto- 
dichroic signal cr • M, where cr is a vector that indicates 
the direction of the incident light and its helicity and M 
is the local magnetization vector. Areas with M paral- 
lel or antiparallel to cr have different secondary electron 
yields, which leads to a bright and dark contrast. In most 
cases, we have chosen an irradiation direction cr perpen- 
dicular to the magnetization direction M in the major 
domains. Although in this case the major domains show 
no magnetic contrast, this choice of the beam direction 
is advantageous because the particularly interesting re- 
gions that separate the domains, i.e. the domain walls, 
display a strong magneto-dichroic signal. The imaging of 
the domain walls requires a very high spatial resolution. 
In the energy range of the Fe 2p3/2 (L3) absortion edge, 
the photon flux in the illuminated area is about 10^^ 
photons per second, which allows for a fast acquisition of 
single images, in the range from 10 s to 30 s for images 
with a field of view of 2.5 fim diameter. The best lat- 
eral resolution in XMCD-PEEM mode has been achieved 
by recording series of about 30 images for each photon 
beam helicity. The images are post-processed using a 
self-correlation algorithm |55j | to minimize drift effects. 
Finally, the magnetic contrast is obtained by subtracting 
graphically the two accumulated images at opposite light 
helicities. The magnetic domain configurations have been 
correlated to the precise morphology of the nanostruc- 
tures using LEEM. The contrast in LEEM is determined 
by the surface topography and its crystalline structure. 
In addition to the diffraction contrast, the so-called in- 
terference contrast allows one to image surface steps and 
thickness gradients in thin films with atomic depth sensi- 
tivity and a lateral resolution of few tens nm. This LEEM 
interference contrast allowed us to unambiguously deter- 
mine the geometry of the self-organized Fe nanostruc- 
tures and of their surrounding 45° inclined facets. The 
quantitative values extracted from the LEEM images (is- 
lands shape, length, width and height) served as input 
for the micromagnetic simulations. A typical LEEM im- 
age and an XMCD-PEEM image of the same sample are 
shown in Fig.|2Ia) and (b), respectively. 



IV. MICROMAGNETICS AND NUMERICAL 
METHOD 

In the framework of micromagnetism the magnetiza- 
tion inside the sample is represented as a directional field 
M = M(r) with the constraint |M| = Ms = const., 
where Ms is the saturation magnetization. For a descrip- 
tion of the intrinsic properties of the magnetic material, 
a set of material constants is used. Besides the satura- 
tion magnetization Mg, the exchange constant A and the 
anisotropy constant K are required. The material param- 
eters are connected with micromagnetic energy terms. In 




FIG. 2: Example of an experimental observation of the par- 
ticle shape and the magnetization structure, (a): Two Fe 
islands, incidentally located next to each other, imaged with 
LEEM. (b): The XMCD-PEEM image displays the in-plane 
magnetization component parallel to the incident beam. The 
beam direction is sketched by the white arrow. Due to grazing 
incidence of the photons (74° with respect to the plane nor- 
mal) the back side of the island is shadowed, thus it appears 
dark in the XMCD-PEEM image. 



our case of magnetic particles with cubic anisotropy in 
absence of an external magnetic field, the relevant energy 
terms are the exchange energy 
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We have used = 4.8 x 10^ J/m^, A = 2.0- 10-"J/m 
and Ms — 1.73 x 10^ A/m to describe the material 
of the Fe islands. The easy axes are oriented along 
the [001], [100] and [010] directions, cf. Fig. [Hd). In 
eq. (O, ai, 012 and as are the directional cosines of 
the magnetization M with respect to the cubic easy 
anisotropy axes. These material parameters lead to a 
quality factor Q = 2K/noM^ ~ 2.5 x 10"^. A value 
of Q ^ 1 indicates that the material is magnetically 
soft, i.e. that the arrangement of the magnetization is 
primarily driven by the need to minimize the der nag ne- 
tizing energy rather than the anisotropy energy. |56| In 
eq. 0, Hs = —'VU is the stray field and U is the mag- 
netic scalar potential. The potential U satisfies Pois- 
son's equation AU = p/^o inside the sample volume V 
and the Laplace equation AU = outside the sample. 
The volume charges p = hq'VM and the surface charges 
a = n ■ M (n: outward surface normal vector) are the 
sources of the stray field. The reduction or avoidance of 
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these charges lowers the stray field energy 0- For the 
inward limit [/; and the outward limit Uo of the potential 
at the surface, the Neumann boundary conditions at the 
surface dV read 



dUj 
dn 



dV 



dUo 
dn 



(4) 



dV 



In addition, the boundary condition lima;^oo U{x) — 
has to be fulfilled. Details on the numerical calculation 
of U by means of a combined boundary element / finite 
element method are given elsewhere p38|. 

If the material parameters and the anisotropy axes are 
specified, the magnetic energy of the sample is uniquely 
determined by the magnetization distribution M{r). In 
the static equilibrium, the magnetization arranges in a 
way to minimize the total energy i?tot = -^cxc + -^cub + 
^'stray In thc computcr simulation, the minimization is 
performed by using the conjugate gradient method [39l| . 
The minimization has to be performed under the con- 
straint \M\ — Ms, which can be easily ensured by repre- 
senting the local magnetization direction with spherical 
coordinates 'd and if. 

The numerical representation and the micromagnetic 
modelling is performed with the finite element method 
(FEM), using a code developed by Hertel ^^J. In thc fi- 
nite element representation, the sample's volume is sub- 
divided into tetrahedral elements of irregular size and 
shape. The magnetization is discretized at the corner 
points of the elements (the nodes). A piecewise lin- 
ear representation of the magnetization is obtained by 
means of normalized, linear interpolation functions that 
are defined in each element (the so-called shape func- 
tions). The FEM is a particularly powerful and accu- 
rate technique for micromagnetic simulations. It allows 
the smooth modelling of particles of general geometry 
(with curved or inclined surfaces) and gives the possibil- 
ity to apply adaptive mesh refinement techniques. These 
unique features of the FEM are of essential importance 
for the present study where large samples with a compli- 
cated shape are modelled. 



V. THE VAN DEN BERG SCHEME 

In thin-film elements, magnetic fiux-closure patterns 
are relatively simple and well-known. Thc Landau pat- 
tern and the diamond pattern (cf. Fig. (SJ are typical 
examples thereof. Such flux-closure patterns are deter- 
mined by the particle's boundary and the tendency to 
avoid magnetic surface charges. In sufficiently large soft- 
magnetic thin-film elements, the magnetization is ori- 
ented in the film plane and is locally aligned parallel to 
thc particle edges. This avoidance of magnetic surface 
charges gives rise to the formation of magnetic domains. 
The Van den Berg scheme |^ ^3 ^ is a simple graph- 
ical method for the construction of flux-closure patterns 
in thin-film elements of arbitrary shape. This method is 
derived from a rather complicated mathematical analysis 




FIG. 3: (color online) Simple magnetic flux-closure patterns 
in thin-film elements, a), c): Landau und diamond structure 
in a rectangle, respectively, b),d): modified Landau and dia- 
mond structure, respectively, in a hexagonal particle. In the 
hexagonal element, the closure domains are subdivided in two 
regions separated by an additional domain wall (dashed line) . 



that assumes an idealized model for two-dimensional thin 
film elements with (a) no magnetic anisotropy, (b) per- 
fect avoidance of magnetic charges, and (c) domain walls 
of vanishing width. According to this scheme, the center 
point of each circle that touches the particle's boundary 
at least twice is placed on a domain wall. This con- 
struction scheme generally leads to symmetric domain 
structures, that, despite the idealizations assumed in the 
model, correspond very well to domain structures as they 
are obtained in experiments 4] and micromagnetic simu- 
lations The Van den Berg scheme was also found 
to be valid in self-assembled Fe/Mo(110) elements of 
moderate thickness (t ^60 nm) |22|. In our thicker sam- 
ples, however, we observe fiux-closure domain structures 
that differ significantly from the Van den Berg scheme. 
It should be noted that the occurrence of deviations from 
the Van den Berg in thick particles is, per se, not surpris- 
ing, since the scheme has been derived using a thin-film 
approximation. A comparable scheme has not yet been 
reported for three-dimensional samples. Analytic studies 
on thc magnetic structure in three-dimensional rnagnetic 
structures have been reported by Arrott et. al. , 
which were, however, restricted to special particle shapes. 
In this paper we report on systematic, qualitative devi- 
ations from the Van den Berg scheme in particles of el- 
evated thickness. In order to highlight these differences, 
we shall refer to the idealized fiux-closure pattern accord- 
ing to Fig|2^) as the classical Landau structure, which 
will serve us as a paradigm for magnetic domain struc- 
tures in thin-film elements. This term is coined with 
reference to the original work by Landau ^45'| , where the 
internal wall structure had been neglected. The classi- 
cal Landau structure is a flux-closure domain pattern in 
rectangular samples. It contains four 90° walls and one 
180° domain wall in the middle. The domain walls are 
assumed to be straight, thin lines. The 90° walls extend 
from the corners to the central 180° wall, which is paral- 
lel to the long edge of the rectangle and with which they 
enclose an angle of 135°. To a very good approximation, 
this pattern corresponds to experimental observations of 
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domain structures in rectangular magnetic platelets. But 
also in particles of non-rectangular shape and even in 
extended films ^^is type of structure has been ob- 
served. In a hexagonal particle, the four-domain Landau 
in a rectangle turns into a structure with six domains 
[cf. FigEb)]. A distinction between this structure with 
six domains and the Landau pattern may be neglected 
since the domain pattern is very similar to the classical 
structure with four domains: two major domains and two 
closure domains. Therefore, such modified variants also 
have been reasonably labeled as "Landau structure" , or 
sometimes as "Landau-type structure" in the literature, 
in spite of more or less pronounced differences from the 
classical structure. In our case of thick particles, how- 
ever, the differences are significant and we require a dif- 
ferent wording. Our study suggests that in the three- 
dimensional case, the classical Landau structure is re- 
placed by a complex flux closure pattern with character- 
istic features. We shall call this the generalized Landau 
structure and we will discuss this structure in detail in 
section IVll Al 



as can be seen by comparing fig. and fig. The 
grey scale in these schematic drawings refers again to the 
y component, so that the dashed domain walls sketched 
in Fig. ^) would not be visible in this experiment, since 
they separate regions of different magnetization direction 
with the same value of my. Similar to the classical Lan- 
dau structure, the magnetization distribution in Fig.^jD) 
is split into four domains: two major domains where the 
magnetization is parallel to the long edge and two closure 
domains that connect the magnetic flux between the ma- 
jor domains. However, neither the shape of the domains 
nor the position of the domain walls corresponds to the 
Van den Berg scheme. 

We have found such distorted, asymmetric domain pat- 
terns bearing only a weak resemblance to the classical 
Landau structure in all the particles that we have inves- 
tigated, as we will discuss in section IVl CI The reason 
for the asymmetric structures and the deviations from 
the Van den Berg scheme lies in the three-dimensional 
nature of these particles. 



VI. RESULTS 

A. Asymmetric flux-closure patterns 

We have taken LEEM and XMCD-PEEM images of 
the magnetization at zero field of about thirty Fe is- 
lands. In all cases, flux-closure magnetization patterns 
have been observed. 




FIG. 4: (color online) Difference between idealized and ob- 
served magnetic domain structure, a) LEEM image of the 
sample shape. XMCD-PEEM image of the magnetic domain 
structure, c) Schematic domain structure for the given shape 
according to the Van den Berg scheme, d) Approximate shape 
of the magnetic domains according to the experimental obser- 
vation. 

Fig.Eb) shows an example of a typical domain pattern 
observed in the experiments. The grey scale displays the 
y component of the magnetization, i.e. the projection of 
the magnetization at the surface along the polarization of 
the incident beam. The magnetic structure in fig. ^p) is 
clearly subdivided into domains, but the shape of the do- 
mains is not in accordance with the Van den Berg scheme. 



B. Micromagnetic Simulations 

The XMCD-PEEM experiments only provide informa- 
tion about one magnetization component (the component 
parallel to the incident beam) on the flat, topmost surface 
of the sample. Therefore, numerical simulation methods 
are employed to investigate the internal structure of the 
magnetization. The simulation results are cross-checked 
with the experimental observations to guarantee that the 
simulated structures correspond to the observed ones. 

Simulating the magnetization in these particles is prob- 
lematic because of their size, which is enormous from the 
point of view of computational micromagnetism. In par- 
ticles of large size, it is difficult to avoid discretization 
errors that result from too large discretization cells. The 
consequence of discretization errors is not only a rough 
representation of the magnetization, but also wrong re- 
sults due to numerical "domain wall collapse" ^3 1 *• ^- 
unrealistic, large change of the magnetization direction 
within a single disretization cell. To avoid discretiza- 
tion errors, the size of the discretiz ation cells sh ould not 
exceed the exchange length A = \/2A/ {poM'^) ~ 3 nm. 
Since the particle volume is of the order of l/im"^, simula- 
tions with a regular grid would require an unsustainably 
large number of several millions of discretization cells. 
An adaptive mesh refinement technique has been 
employed to solve this problem. Such adaptive methods 
are particularly powerful for static micromagnetic sim- 
ulations of large particles, where the inhomogeneities of 
the magnetization are limited to a small fraction of the 
sample's volume. These regions, i.e. the domain walls, 
require a high discretization density, whereas large do- 
mains with homogeneous magnetization can be calcu- 
lated safely with a coarse mesh. Therefore, adaptive re- 
finement methods allow for accurate numerical results 
while keeping the computational costs (processing time. 
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FIG. 5: (color online) To avoid discretization errors in the 
simulation of the three-dimensional magnetic domain struc- 
ture, the sample is first reduced to a tenth of its size. By 
using cycles of micromagnetic energy minimization and three- 
dimensional adaptive mesh refinement, the sample is gradu- 
ally increased to its real size (a). After this procedure, a finite- 
element mesh with regions of strongly different discretization 
density is obtained (b). The finite elements are particularly 
small in regions of strong magnetic inhomogeneities (vortices, 
domain walls). Owing to the mesh refinement, even the nano- 
metric vortex core with magnetization perpendicular to the 
surface is fully resolved, as can be seen in the red spot in 
panel (c). Throughout the paper, unless otherwise specified, 
the color coding used to represent a magnetization component 
is chosen according to the bar on the bottom of panel (a). 

memory) low. 

The procedure we have applied for the micromagnetic 
simulations is as follows. First, the exact geometrical 
shape of each Fe island is extracted from LEEM images 
and a corresponding finite element model is constructed. 
Subsequently, the model is scaled down to 1/10 of its ac- 
tual size. The magnetization distribution in this small is- 
land is simulated by means of energy minimization, using 
a symmetric vortex structure around the center as initial 
configuration. Once a converged solution is found, the 
mesh is adaptively refined and the equilibrium structure 
is calculated with the refined mesh. The cycle of mesh 
refinement and energy minimization is repeated until the 
maximum angle enclosed between the magnetization vec- 
tors at two neighboring discretization points drops below 
7r/10. Then, the size of the sample is increased, and the 
procedure of mesh refinement and energy minimization 
is repeated. With this iterative process, the sample size 
is increased in small steps towards its real value and the 
mesh is progressively refined. This way of slowly "in- 
flating" the sample guarantees that the magnetic struc- 
ture is calculated without discretization errors. The sim- 



pler variant of simulating the magnetic structure directly 
with a coarse mesh, without reducing the sample's size, 
and subsequently applying the mesh refinement proce- 
dure is less accurate because the first calculation would 
lead to large discretization errors. Especially in three- 
dimensional calculations, it is advisable to prevent dis- 
cretization errors rather than trying to remove them a 
posteriori. An example for the simulation of the mag- 
netization by means of stepwise increase of the sample 
size is shown in Fig. ]^) . The typical cell size of a re- 
fined mesh is between about 25 nm in the homogeneous 
regions and 1.5 nm in the refined regions. By increasing 
the sample size, we find that the symmetric vortex struc- 
ture in small particles evolves to an asymmetric domain 
structure above a certain sample size. Such a transition 
can be seen in Fig. where the asymmetry occurs af- 
ter the sample has been increased to 50% of its actual 
size. The formation of the asymmetric structure obvi- 
ously involves a breaking of symmetry: the vortex which 
is originally placed in the center may either shift to the 
right or to the left side on the top surface. As will be 
discussed later, the sign of the perpendicular component 
in the vortex core may also have a strong influence on 
the resulting domain structure. 

C. Comparison between Experiment and 
Simulation 

The simulation results are in excellent agreement with 
the experimental data. A number of examples are shown 
in Fig.|Hl The experimentally observed asymmetric shape 
of the end domains and the position of the domain walls 
are reproduced almost perfectly by the simulations. All 
the unexpected features of the domain patterns and their 
deviation from the Van den Berg scheme occur in just the 
same way in both the experiment and the simulations. 
Notice that all these details of the magnetic structure 
develop automatically in the simulation by using a simple 
vortex as a starting configuration. 

The one-to-one correspondence between experimental 
and simulated surface magnetization is a strong indica- 
tion that the simulated structures are indeed equal to 
the experimental ones. Therefore, the simulations can be 
used to obtain information that is not directly accessible 
in the experiment, i.e. the three-dimensional structure of 
the magnetization inside the sample. The investigation 
of the inner structure of the magnetization is a prerequi- 
site for the explanation of the unexpected experimental 
results. 



VII. DISCUSSION 

A. Internal domain wall strucure 

The magnetization distribution inside the sample dif- 
fers significantly from the distribution observed at the 
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FIG. 6: Experimental and simulated data for five self-assembled Fe/Mo(110) islands of different size and shape (a-e). First row: 
The shape of the particle is obtained by LEEM imaging. The particle thickness is derived from the width of the (001) facets, 
which appear as a dark border around the hexagonal top surface, by using the known inclination angles of the inclined facets 
(cf. Fig. 0. The sample's approximate lateral dimension and thickness is displayed on top of each column. Second row: One 
magnetization component in the surface plane is imaged as grey scale with XMCD-PEEM. In sample (a), the contrast refers 
to the magnetization component parallel to the long edge while in the other samples the in-plane magnetization component 
perpendicular to the long edge is displayed. Notice that only the top surface, i.e. the light grey internal hexagon in the first 
row, is imaged. Third row: Micromagnetic simulation results. To compare the results with the XMCD-PEEM experiments, 
the same magnetization component is displayed in grey scale. The sample shape was modelled according to the LEEM images 
and the model explained in Fig. 



surface. A typical example is shown in Fig.[7^),b,)c) that 
displays the y component of the magnetization on three 
cross-sections at the top {z = t), the middle {z = t/2) 
and the bottom (z — 0), respectively. The coordinate 
frame is chosen according to Fig. ^ In the middle cross 
section, the domain pattern is almost symmetric and in 
good agreement with the Van den Berg scheme, while on 
the top and the bottom the structure is clearly asym- 
metric. The twist is in opposite directions on the top 
and the bottom surface, cf.[7^),c). The inhomogeneity of 
the magnetization along the z-direction and the resulting 
asymmetry in the top and bottom xy-plane result from 
the formation of an asymmetric 180° Bloch wall [26ll4^ : 
in the middle of the sample (z — t/2) the major domains 
are separated by a Bloch wall, while on the top and bot- 
tom surface, the transition between the major domains is 
given by a Neel wall (the so-called Neel caps) . In a cross 
section on a plane with x =const. through the middle 
of the sample, the typical profile of an asymmetric Bloch 
wall resulting from the combination of Bloch and Neel 
wall can be clearly seen[7|i.),i),j). There are different pos- 
sibilities for the orientation of the magnetization in this 
type of domain wall: The magnetization in the Bloch 
part can point either in the positive or in the negative 
z direction, and the magnetization in the Neel caps can 
point either in the positive or negative y direction. The 
orientation or the chirality of the asymmetric Bloch wall 
can have a strong impact on the resulting domain struc- 
ture, as shown in Fig. |S1 Depending on the direction of 
the magnetization of the Neel caps and the direction of 
the inner Bloch component, the magnetic domain struc- 
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FIG. 7: (color online) Magnetization distribution on differ- 
ent cross sections. The example shown here refers to the 
particle displayed in Fig. (SJj) . The figures (a) and (d) show 
the top view on the three-dimensional sample with lateral 
facets. (b),(e): cross-section through the middle of the sam- 
ple at z = t/2; (c),(f): Magnetization on the bottom surface 
of the sample, z = 0. (g): component on the cross-section 
through the middle of the sample parallel to the xz-plane. 
(h-j): asymmetric Bloch wall: magnetization components on 
a cross-section at mid-length, parallel to the yz-plane. 



ture on the surface may be very different. In the example 
shown in Fig.|Sl the computer simulation first yielded the 
structure (a) for this geometry. In the simulations, the 
aforementioned breaking of symmetry that is involved 
with the formation of such structures results from numer- 





FIG. 8: (color online) Simulation of four different configura- 
tions in a sample of identical shape. The shape and the size of 
the sample is the same as in fig. |SJi) . Although topologically 
the structure is very similar in all cases (Landau structure 
with asymmetric Bloch wall) , the resulting magnetic patterns 
are significantly different. In the structures (a) and (b) the 
Bloch component of the internal wall, m^, is positive, whereas 
in (c) and (d) it is negative. The structures (a) and (c) have 
the same direction of the magnetization in the Neel caps, my, 
which is opposite to the orientation in the structures (b) and 
(d). 



ical roundoff errors and the direction in which the Neel 
caps evolve is impredictable. This imponderability can 
be removed by breaking the symmetry with an oblique, 
weak external field that is applied during the simulated 
expansion of the sample, cf. Fig. (SJi). With this con- 
trolled breaking of symmetry we obtained the structures 
shown in Fig. |H1 including the structure (d), which al- 
most perfectly reproduces the experimental observation, 
cf . Fig. IHli) . Apart from possible minor differences con- 
cerning the out-of plane component of the magnetization 
along the facets, there are four main types of such almost 
degenerate structures, corresponding to the two possibil- 
ities for the orientation of the Bloch component and the 
Neel caps, respectively [s^- In a sense, these structures 
can be compared to geometric isomers of a molecule: The 
domain structures contain the same components (Lan- 
dau structure with asymmetric internal walls), but the 
orientation is different (opposite orientation of the Neel 
caps and/or of the Bloch wall). The domain structures 
of such "isomers" can differ significantly in their patterns 
on the surface, especially if the sample shape is strongly 
asymmetric. 

In any case, the Neel caps of an asymetric Bloch wall 
are oriented in opposite direction on opposite surfaces. 
The formation of an asymmetric Bloch wall represents an 
energetically optimized arrangement concerning the stray 
field energy. Magnetostatic volume charges p = V • M 
are largely avoided by the formation of a Bloch wall in 
the center of the sample and surface charges a = M ■ h 
(fi: surface normal vector) are effectively suppressed by 
the Neel caps. Depending on whether the magnetization 
in the Neel caps is pointing parallel or antiparallel to 
the magnetization direction in the closure domains (the 



FIG. 9: (color online) Geometry dependence of the magne- 
tization structure in a three-dimensional island, (a): Per- 
spective view on the simulated magnetization structure of the 
particle shown in Fig. |S|d). (b): Simulated magnetic domain 
structure in a particle of same size and thichkness, but with 
symmetric in-plane shape and vertical facets, (c): Even if the 
hexagonal shape is replaced by a smooth elliptical boundary, 
the magnetization structures remains almost unchanged, with 
the same asymmetric characteristics. 



domains close to the particle's border in the ±a; direc- 
tion), the Neel caps are either connected smoothly to the 
closure domain or separated by them with a magnetic 
vortex, respectively 50]. Inside the closure domain, on 
cross-sections with constant value of x, the structure is 
similar to an asymmetric Neel wall. The main differ- 
ence between an asymmetric Neel wall and an asymmet- 
ric Bloch wall is the relative orientation of the magne- 
tization in the Neel caps, being parallel in the former 
and antiparallel in the latter. It has been predicted, by 
means of numerical simulations, that the combination of 
asymmetric wall types in mesoscopic, soft-magnetic par- 
ticles leads to a pro nounced in-plane asymmetry on the 
top surfaces HI 113. Our XMCD-PEEM investigations 
on several mesoscopic Fe islands clearly confirm this pre- 
diction and provide the first direct observation of this 
effect. 



B. Shape dependence 

In the numerous particles that we have investigated, 
we have almost exclusively found the Landau-type flux- 
closure pattern described above, the generalized Landau 
structure. A few exceptions of different magnetization 
arrangements will be discussed in the next sections. The 
generalized Landau structure consists of four domains, 
two major domains that are connected by an asymmet- 
ric Bloch wall and two closure domains connected to the 
major domains by 90° walls. Our observation of the Lan- 
dau structure in self- assembled Fe islands with their quite 
complicated hexagonal shape and the inclined facets sug- 
gests that the formation of this arrangement of the mag- 
netization is rather driven by the particle's aspect ratio 
and its size than by the rectangular shape. 

The question arises to which extent the asymmetric 
Landau pattern that we observed is specific to our ex- 
ample of hexagonal particles with inclined facets. How 
sensitively does the magnetic structure depend on the 
precise shape of the particle? To answer this question, 
the influence of the particle shape on the resulting mag- 
netic domain structure has been investigated with micro- 
magnetic simulations. Unlike the real world, where the 
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shape of the Fe islands is determined by the energetics 
of the growth mode, finite element models give the pos- 
sibility to study magnetic structures in mesoscopic par- 
ticles of arbitrary shape. Starting from a "real" island, 
i.e. one that is modelled exactly according to the exper- 
iment, we have first tested to which extent the magnetic 
pattern is affected by the inclined facets and the asym- 
metry of the hexagon. This has been done by simulating 
the magnetization in a hexagonal particle of same thick- 
ness and size, but with symmetric hexagonal shape and 
a "cookie-cutter" geometry, i.e. an island with perpen- 
dicular facets. The result is displayed in Fig. O)). The 
asymmetric magnetic structure is practically the same 
as it has been observed in the real island with inclined 
facets, indicating that neither the inclination of the facets 
nor the asymmetry of the hexagon has a decisive impact 
on the magnetic structure. To further check whether the 
hexagonal shape is of importance for the resulting mag- 
netic structure, a flat, elliptical particle has been mod- 
elled. The major axes of the ellipse (parallel to the x 
and y direction, respectively) correspond to the width 
and length of the original island and the thickness is 
again the same. The result [Fig. ^)] shows that the dis- 
torted Landau structure also occurs in this particle. Re- 
cent studies show that practically the same structure 
may also be found in an ellipsoidal particle of compara- 
ble size and aspect ratio. These computer experiments 
indicates that the generalized Landau structure is hardly 
determined by the particle shape. This is in contrast 
to two-dimensional samples, where the essential impor- 
tance of the sample shape reflects in the Van den Berg 
construction scheme, where the domain structure is de- 
termined exclusively by means of the particle shape. The 
occurrence of the three-dimensional, generalized Landau 
structure is neither restricted to particles of rectangu- 
lar, hexagonal or any other speciflc shape; not even to 
flat particles. In the three-dimensional case, it is rather 
the thickness, the size and the aspect ratio that are de- 
cisive. The primary tendency of the magnetization in 
thick, three-dimensional particles appears to be to de- 
velop an asymmetric 180° Bloch wall, while in the case 
of thin-film elements, the major domains may either be 
separated by a 180° wall (classical Landau structure) or 
by 90° walls (diamond state). 



C. Diamond structure 

By far, the magnetization states that we have observed 
most frequently are the above-mentioned variants of the 
Landau structure. The absence of the "diamond state" 
in these thick elements is remarkable. The diamond state 
is an alternative flux-closure pattern with seven domains 
(cf. FigO that is frequently observed in thinner soft- 
magnetic particles [23,|33,|51| of similar size. 

We suspect that the predominance of the Landau 
structure and the absence of the diamond state can be ex- 
plained with energetic considerations in connection with 




FIG. 10: In a relatively large and particularly thin island, 
identifiable as such in a LEEM image (a) by the thin dark 
boundary (the lateral facets), a "diamond state" magnetiza- 
tion structure was found. The diamond state is character- 
ized by two vortices with opposite vorticity and a diamond- 
shaped central domain, which is clearly recognizable on the 
right panel (b). 

the elevated particle thickness. The thickness depen- 
dence of domain configurations in nanoscale Fe/W(110) 
islands of significantly lower thickness (ca. 3-8 nm) has 
been studied recently by Bode et al. . In the present 
case, the idea is the following. Essentially, the Landau 
structure has one central 180° wall and four 90° domain 
walls, whereas the diamond structure consists of seven 
domains that are separated by eight 90° domain walls. 
While in the case of a 180° wall magnetic volume charges 
can efficiently be avoided by the formation of an asym- 
metric Bloch wall, a 90° wall inevitably contains volume 
charges p = VAf, like a Neel wall does. With increas- 
ing thickness, the energy connected with these volume 
charges increases and so does the tendency to avoid them. 
The tendency to avoid these charged 90° walls and, in- 
stead, to form a 180° wall, favors the formation of the 
Landau structure over the diamond structure. If this 
consideration is correct, the tendency to avoid the di- 
amond state should decrease with decreasing thickness. 
We therefore searched for the diamond state in islands 
of large in-plane extension (to ensure that the particle is 
in a multi-domain state) that were particularly thin. In 
our set of samples the typical thickness was about 100- 
150 nm, but we eventually found an island of 65 nm thick- 
ness and about 1000 nm x 500 nm lateral extension. The 
XMCD-PEEM investigation of the magnetization in this 
island indeed showed that the sample was magnetized in 
the diamond state, cf. Fig. ^| This sample is the only 
one in which we found the diamond state. The fact that 
we have observed this state only after a selective search 
aimed at finding particularly thin elements suggests that 
the thickness, actually, plays a decisive role. Further 
studies are required to systematically analyze this sug- 
gested thickness dependence of magnetization states. 



D. Vortex structure 

Besides the thickness and the size, also the in-plane 
aspect ratio has an impact on the resulting magnetiza- 
tion state. This can be seen by comparing the magnetic 
structure of two islands of similar thickness and width, 



but diflFerent length, cf. Fig. |Hl)b),c). In one case, the 
aspect ratio (length in x direction vs. width in y direc- 
tion) is 2.1 and in the other it is about 1.5. The particle 
with lower aspect ratio has a vortex in the center of the 
sample. It contains no asymmetric Bloch wall and the 
closure domains are asymmetric and twisted. They are 
stretched in the x direction and they join in the middle of 
the sample at the vortex. In contrast to this, the longer 
sample has a clearly developed 180° domain wall in the 
middle which separates the major domains. Obviously, a 
sufhciently large particle with aspect ratio approximately 
equal to one would not form a Landau structure with a 
central 180° domain wall, but would be magnetized in a 
simple vortex state with a magnetic vortex structure in 
the centre of the sample. With increasing aspect ratio, 
the vortex core is stretched until it eventually performs a 
kink through the sample and forms the center of an asym- 
metric Bloch wall. Thus, by changing the aspect ratio, 
a transition from a vortex state to a generalized Landau 
structure occurs at a critical size. Whether this transi- 
tion is continuous or discontinuous may be the subject of 
future studies. 



E. Bloch switch 

In most cases the simulations yield magnetic structures 
that agree very well with the ones that have been ob- 
served in the experiments, but occasionally some interest- 
ing differences between simulation and experiments have 
occurred. The possible differences concerning the chi- 
rality of the asymmetric central Bloch wall have already 
been discussed in section IVII Al Another variant of the 
Landau structure was obtained in one case, where the 
simulation has yielded a structure as shown in Fig. lllb ). 
c). Instead of the usual straight line, the central do- 
main wall (i.e. the Neel cap) displays a pronounced kink. 
It is known from studies of asymmetric domain walls in 
extended magnetic films that such a kink is characteris- 
tic for a so-called Bloch switch |^ [s^ . A Bloch switch 
is given when the sign of the Bloch component changes 
along the domain wall inside the film, while the magneti- 
zation direction of the Neel caps remains unchanged. In 
fact, a look at the inner part of the simulated magneti- 
zation structure [Fig. IllF)] shows that the arrangement 
contains Bloch switches. We have also observed a similar 
kink in the experiment, as shown in Fig. 1111 ;). 

This example is the only island in which we have found 
this kink experimentally, and it can therefore not be ruled 
out that the kink may have other reasons than a Bloch 
switch like, e.^., a structural inhomogeneity of the island. 
Moreover, the top surface of the island in which we ob- 
served the kink was not ideally flat, as could be seen by 
the shadow of the island as described in Fig.j^ti). Due to 
these deviations of the sample morphology from the ideal 
construction scheme according to Fig. ^ the model for 
the island shape that has been used in the simulations 
that yielded a Bloch switch does not correspond exactly 
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FIG. 11; (color online) (a-c): z,x and y component of the 
magnetization of a simulated magnetic domain structure with 
undulated central domain wall, (f) ^-component of the mag- 
netization on a cross-section of constant y. The alternating 
red- blue regions indicate switches of the Bloch component, 
(d): LEEM image and XMCD-PEEM image (e) of a sample 
in an apparently similar magnetization state. The size of the 
island is 1500 nm x 620 nm x 130 nm. 



to the island in which the kink has been observed experi- 
mentally. Therefore, although a correspondence between 
experiment and simulation is likely, the available data 
from simulation and experiment are not sufficient for an 
unambiguous proof that the kink results from a Bloch 
switch. The results, however, suggest that an asymmet- 
ric Bloch wall with a Bloch switch may develop also in 
patterned particles and that such a Bloch switch would 
lead to a considerable distortion of the otherwise straight 
central domain wall. To our knowledge, this kind of do- 
main wall has not been reported previously in patterned 
elements. 



VIII. CONCLUSIONS 

"The simplest problem in continuum magnetism in 
a singly connected finite body is the Landau structure 
(...), but it remains an unsolved problem to describe 
it in all its details." This is a quote from a paper by 
A.S. Arrott and T.L. Templeton ^O] that was written 
in 1997, at a time when the experimental resolution was 
not yet sufficient for a direct observation of details of the 
magnetization and most micromagnetic simulation meth- 
ods were not accurate enough for large, three-dimensional 
computations |59j . Owing to the combination of (a) ad- 
vanced growth methods yielding high-quality, monocrys- 
talline mesoscopic Fe particles that can be regarded as 
model systems, (b) LEEM imaging to obtain detailed 
information on the particle shape, (c) high-resolution 
XMCD-PEEM studies that provide information on the 
in-plane manetization at the particle surface, and (d) ac- 
curate FEM/BEM micromagnetic modelling with adap- 
tive mesh refinement techniques providing a detailed 
model of the three-dimensional magnetization structure, 
we have solved the problem of a full description of the 
Landau structure. The Landau structure seems to occur 
as a minimum energy arrangement in soft-magnetic parti- 
cles of very different shape. Generally, the size, thickness 
and aspect ratio are more important for the formation of 
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the generalized Landau structure than details of the par- 
ticle's shape and its surfaces. The fascinating complexity 
of the Landau structure occurs when the particle is thick 
enough to sustain an asymmetric Bloch wall. The classi- 
cal, symmetric Landau pattern occurring in rectangular 
thin film elements can be regarded as a simple thin-film 
limit of a more general, three-dimensional structure. In 
the thin- film limit, the central asymmetric Bloch wall 
is replaced by an ordinary Neel walls, thus simplifying 
drastically the arrangement. 
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